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Information, disturbance and Hamiltonian quantum feedback control 
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We consider separating the problem of designing Hamiltonian quantum feedback control algo- 
rithms into a measurement (estimation) strategy and a feedback (control) strategy, and consider 
optimizing desirable properties of each under the minimal constraint that the available strength 
of both is limited. This motivates concepts of information extraction and disturbance which are 
distinct from those usually considered in quantum information theory. Using these concepts we 
identify an information trade-off in quantum feedback control. 
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I. INTRODUCTION 



With experimental advances, particularly in the fields 
of cavity QED § and ion trapping [Q, it is possible to 
observe individual quantum systems in real time, and it 
is therefore natural to consider the possibility of control- 
ling such systems in real time using feedback [pj. As a 
special case, real-time Markovian quantum feedback has 
been analyzed (J |7j and implemented experimentally in 
certain quantum optical systems ||. Feedback control 
is invaluable in macroscopic applications, and as a con- 
sequence there is a vast body of literature devoted to 
classical feedback control. While results from classical 
control theory may be applied fruitfully to the quantum 
domain in certain special cases |||l0), these are not ad- 
equate in general, primarily because quantum measure- 
ment is quite different in nature from classical measure- 
ment, in that it has the capacity to disturb the system 
under observation PJ. As a result, the development of 
optimal quantum control strategies requires optimizing 
over possible measurement strategies, which is unneces- 
sary in classical control. 

In feedback control the dynamics of a system is manip- 
ulated by using information obtained about the system 
through measurement. The goal is usually to maintain 
a desired state or dynamics in the presence of noise. A 
central problem of feedback control theory is the devel- 
opment of algorithms to achieve this goal. The approach 
to controller design that we consider here is to examine 
the measurement and feedback steps separately, thereby 
splitting the feedback control problem into two parts. 
One can then consider optimizing desirable properties of 
these parts separately under suitable constraints. If one 
allows the strength of either measurement or Hamilto- 
nian feedback to be infinite, then any control objective 
can be achieved perfectly (this will be shown below once 
we have made these concepts of strength more precise) . A 
constraint on strength is therefore the minimal constraint 
under which the problem of quantum feedback control is 
non-trivial, and this is the constraint we employ here. 

The action of optimizing for the feedback and mea- 
surement independently ignores the possibility that truly 



optimal solutions may require considering both together. 
We will also simplify the problem by considering the op- 
timization at each time step separately. This assumes 
that it is never desirable to perform worse at the cur- 
rent time in order to perform better at some future time. 
The approach we take here is therefore not aimed at find- 
ing a globally optimal solution given a set of constraints. 
However, the expectation is that the concepts we intro- 
duce here provide a simple systematic approach which 
one can expect to produce good results, and provide an 
insight into the kind of measurement processes which are 
desirable in feedback control. 

For the feedback step, we consider the question of the 
effectiveness of the control by defining a cost function. 
Since one is interested in controlling the dynamics of a 
given quantum system (usually in the presence of some 
unavoidable source of environmental noise), one can spec- 
ify the objective by specifying the most desired state for 
the system at each instant. The 'cost' function is then 
the sum of the distances of the state of the system from 
the desired state at each point in time, for some suit- 
able measure of distance. Finding the optimal control 
strategy then consists in minimizing the cost function, 
under suitable constraints for the strength of the feed- 
back. Note that this is somewhat different from the 
standard approach taken in modern classical control the- 
ory fl3|-|l5|| , and more similar to the approach taken in the 
new techniques of 'postmodern' classical control In 
modern classical control (e.g. LQG theory) one usually 
optimizes a 'total' cost function obtained from a suitably 
weighted sum of the cost function defined here, and an- 
other cost function intended to capture the cost of feed- 
back strength. 

We will restrict ourselves to control objectives such 
that the desired state at each time (the target state) is 
pure. Impurity (mixing) merely signifies that our knowl- 
edge of the quantum system is less than maximal, which 
is by assumption undesirable. 

In considering the optimality of the measurement step, 
rather than attempting to find a measurement which ex- 
plicitly optimizes the cost function, we define concepts of 
information and disturbance, motivated by the feedback 
control problem. We then consider finding measurements 
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which maximize the information and minimize the dis- 
turbance. We find that in general these two targets are 
mutually exclusive, in striking contrast to classical con- 
trol theory. This implies the existence of a trade-off be- 
tween information and disturbance in quantum feedback 
control. 

Since we focus on continuous feedback control, and 
many readers will be familiar with generalized measure- 
ments but unfamiliar with the formalism of continuous 
quantum measurement, we describe in the next section 
how continuous observation is formulated within the lan- 
guage of generalized measurements. In Section [II we 
define the concept of the strength of a measurement, re- 
quired as a minimal constraint for the feedback control 
problem. In Section |v| we discuss in detail the division of 
feedback control into 'pure' measurement and Hamilto- 
nian feedback, and consider what may be achieved when 
there is no limitation on the strength of either. We also 
discuss what may be achieved in this case both with- 
out feedback and with measurement-only feedback. In 
Section [y], we consider the measurement process, define 
concepts of information and disturbance, and consider 
minimizing the disturbance and maximizing the informa- 
tion. In Section VI we examine the Hamiltonian feedback 



!^ -(i)" 4 - wa -'" , ' 



(i) 



and obtain Hamiltonians which min imize the instanta- 
neous cost function. In Section VII we implement the 
feedback control of a two-state system, showing how the 
ideas presented in the previous sections mani fest i n the 
performance of the control algorithm. Section VIII sum- 
marizes and concludes. 



II. CONTINUOUS OBSERVATION AND 
GENERALIZED MEASUREMENTS 

We will concern ourselves primarily with continuous- 
time quantum feedback control, in which a system is ob- 
served continuously, and the results of the measurements 
(the measurement record) used to continuously alter the 
Hamiltonian of the system to effect control. We now 
discuss how continuous observation may be described 
within the language of generalized quantum measure- 
ments, implemented as positive operator valued measures 
(POVM's). 

Continuous measurements on a quantum system gen- 
erate a measurement record that is a continuous-time 
stochastic process, which may be either a (Gaussian) 
Wiener process or a Poisson process For a given 

physical system, these two kinds of processes will result 
from making different measurements, for example photon 
counting and homodyne detection of optical beams. 

The key ingredient in describing continuous measure- 
ments is that during an infinitesimal time step dt, the 
information obtained by the observer must scale as dt, so 
that one can take the continuum limit and obtain a sen- 
sible answer. This may be realized by defining a POVM, 
given by f W a Q a da = 1, to describe the result of an ob- 
servation in the time interval dt by [H8| 



where Q is a arbitrary operator for the system under ob- 
servation, a takes all values on the real line, and k is 
a positive real constant. For reasons that will be made 
clear in the next section, we will only need to be con- 
cerned with the case in which Q is Hermitian, so that Q 
may be referred to as an observable, and we will assume 
this in what follows. Note that each tt a is a weighted sum 
of projectors onto the eigenbasis of Q, where the weight- 
ing is peaked at a. Thus each application of the one of 
the £l's provides some information about the observable 
Q. However, as dt tends to zero, this information also 
tends to zero, since the fi's become increasingly broad 
over the eigenstates of Q. Calculating the measurement 
result in the interval dt at time t, and denoting this as 
dy(t), we have 



(2) 



dy(t) =4k{Q)dt + V2kdW, 



where dW is the Wiener increment for the interval dt. 
Using this, one can obtain the stochastic evolution of 
the quantum state under this measurement process, re- 
ferred to as a quantum trajectory, and this is given by 
the Stochastic Master Equation (SME) 



dp = -i[H, p]dt - k[Q, [Q, p]]dt 

+ (Qp + PQ- 2Ti[Qp]p)V2kdW. 



(3) 



where H gives the system evolution in the absence of 
the measurement. We can also readily obtain the non- 
selective evolution, in which the measurement results are 
ignored, and this is given by 



p(t + dt) = -i[H, p]dt + J Qapfl^da 

= -i[H,p]dt-k[Q,[Q,p{t))]dt. 



(4) 



When H commutes with Q this evolution leads to a diag- 
onalization of p in the basis of Q, as one would expect for 
a measurement of Q. Similarly, integrating the SME in 
this case, one finds that the result in the long time limit 
is a projection onto one of the eigenstates of Q. Such a 
POVM realizes a continuous measurement of the oper- 
ator Q, such that the measurement record is a Wiener 
process. 

One can also define a POVM to provide continuous ob- 
servation in which the measurement record is a Poisson 
process. Since this requires only one of the two possible 
outcomes at each interval dt, the POVM consists of only 
two measurement operators: 



fi = 1 - \kQ 2 dt, 
Vt x = Q\fkdi. 



(5) 
(6) 



That this gives a Poisson process can be seen by consid- 
ering the probabilities for the outcomes and 1, which 
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are 1 — k(Q)dt and k(Q)dt, respectively. Result 1 there- 
fore corresponds to an Poisson 'event', which happens 
occasionally, and to the absence of one. The SME cor- 
responding to the measurement process is different from 
that corresponding to the Wiener measurement, but the 
non-selective evolution is identical. Physically, the non- 
selective evolution is fixed by choosing the interaction 
of the system with the environment which is mediating 
the measurement, and the trajectory, whether Poisson or 
Wiener, is selected by how one chooses to measure the 
environment so as to extract the information about the 
system. In fact, by taking a suitable unitary transforma- 
tion of the Poisson measurement operators, and taking 
the appropriate limit in which there are many events in 
each interval dt, one can obtain the Wiener process mea- 
surement from the Poisson measurement, and so the first 
can be regarded as a special case of the second. This is 
discussed in detail in |l2|,|l9| . 

The point we wish to note here is that regardless of 
how one chooses the trajectory, a continuous measure- 
ment of an observable Q is given by a POVM in which 
all the measurement operators Q a are positive operators, 
diagonal in the basis of Q, and one must merely be care- 
ful to choose the form of these operators with respect to 
dt so as to provide a sensible continuum limit. 



III. THE STRENGTH OF A MEASUREMENT 

Clearly the more accurate the measurements of the 
observer, the more information she is able to obtain, and 
better able she is to choose feedback to effectively con- 
trol the system. However, in general, more accurate mea- 
surements require more resources. In treating quantum 
feedback control, it is sensible to consider a restriction on 
available resources, and hence a restriction on measure- 
ment accuracy. To treat this quantitatively, one must 
introduce a sufficiently precise notion of the accuracy, or 
strength, of a quantum measurement. 

For the purposes of feedback control, since it is the fi- 
nal state resulting from measurement that the observer 
must act upon with feedback, it is the observers infor- 
mation about this final state which is relevant. Intu- 
itively, one can therefore think of stronger measurements 
as providing, on average, final states which are more pure 
(or, alternatively, have a smaller von Neumann entropy) 
than weaker measurements. When considering continu- 
ous observation, in the absence of any noise sources, an 
initially impure state is continually purified. In this case 
the strength of the measurement can be thought of as be- 
ing proportional to the rate of this purification. Note that 
this concept of information extraction by a measurement 
is quite different from that usually considered in quan- 
tum information theory. There, authors have been con- 
cerned about the information that a measurement pro- 
vides about the initial state of the system (the state im- 
mediately before the measurement) pqjiil{| , whereas in 



our case it is the information about the final state which 
is important. 

We will not need an explicit definition for measurement 
strength here, since we will only require two properties 
of measurement strength which we will motivate below. 
However, we will give an example of an explicit definition 
which satisfies these two properties. To motivate the first 
property, we note that as we have defined it so far, it is 
clear that the strength of a measurement in some sense 
characterizes the average rank of the operators f2 m which 
make up the associated POVM (J2 m ^ti^m = !)• If all 
the £l m are rank one, then one always obtains a pure final 
state, and therefore complete information, regardless of 
the initial state. The higher the rank of the projectors, 
the higher in general will be the von Neumann entropy 
for a fixed initial state. The first property we will require 
is that measurements that consist of rank one projectors 
should have maximum strength (for measurements on a 
system of a given dimension). 

For the remainder of this paper, we will refer to mea- 
surements for which at least one of the f2 m are rank one 
as infinite strength measurements. This terminology is 
natural in the context of continuous observation, since 
in order to provide rank one projections in a finite time 
from a continuous measurement, one would have to take 
the limit k — > oo in Eq.(^). However, we wish to stress 
that our use of this terminology is not intended to im- 
ply that any explicit definition of measurement strength 
should necessarily take this value for these kinds of mea- 
surements. 

The second property we wish to impose is that strength 
be invariant under unitary transformations of the mea- 
surement operators. To motivate this property, one can 
consider a device which measures the spin of a two state 
system. One would expect such a device to provide the 
same strength of measurement regardless of how it is ori- 
ented in space. Since spatial rotation covers all unitary 
transformations for a spin-half, for this system strength 
should be invariant under all unitary transformations of 
the r2„. We will explicitly consider the spin- half system 
later. 

To provide an example of an explicit definition of mea- 
surement strength for single-shot measurements on finite 
dimensional systems, one can first consider the average 
uncertainty after the measurement result is known. Us- 
ing the von Neumann entropy, for a measurement de- 
scribed by J2 n ^n^n = 1, this is 



(p) = T^npnjMfinpfiJ, /Tr[n nP ni ])], (7) 



where p is the initial state of the system. Using the purity 
as an alternative measure of uncertainty we have 



*p(p) = i-E 



Tr[(0^t) 2 



Tjr[O n pfiJ,] 



ti ■ 



(8) 



To obtain definitions of measurement strength that sat- 
isfy our two properties we can use the following functions 
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of these uncertainties: 



1 



1 



U V (I/N) 
1 



log(iV)' 
N 



u p (I/N) (N-iy 



(9) 
(10) 



in which N is the dimension of the system being mea- 
sured. These particular definitions have the additional 
property that they tend to infinity for von Neumann mea- 
surements, but no others. It is also a fairly simple matter 
to write explicit definitions of measurement strength that 
tend to infinity for measurements that contain at least 
one rank one projector. 

Definitions of measurement strength for single-shot 
measurements may be extended to continuous measure- 
ments by using the rate of uncertainty reduction. Using 
the explicit definitions given above (Eq.(g) and Eq^)), 
it is straightforward to calculate this rate from Eq.(m and 
the Ito rules for stochastic differential equations |22]: 
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d 



TVlog(iV) 2 
'(N-l) 



(Tr[Q 2 ]+3Tr[Q] 2 



(11) 
(12) 



In considering feedback control, measurement strength 
is a particularly important concept because it is a con- 
strained resource; stronger measurements are in general 
more expensive. A particular example is the measure- 
ment of position by the reflection of a laser-beam |l0|,^3| , 
a technique used in the atomic force microscope. In 
that case it is the laser power on which the measurement 
strength depends. Hence it is reasonable to consider opti- 
mization of the measurement under the assumption that 
measurement strength is fixed. 



IV. MEASUREMENT AND FEEDBACK 

In classical feedback control, it is natural to consider 
the measurement process as being qualitatively different 
from the feedback process. In particular they may be dis- 
tinguished by the fact that the measurement in each time 
step involves no change to the system Hamiltonian, and 
the feedback step provides no information. In quantum 
feedback, since measurement has the ability to affect the 
dynamics in ways that in classical mechanics would have 
to be attributed to a Hamiltonian, the distinction is not 
as fundamental. However, in the vast majority of quan- 
tum feedback schemes considered to date, it is some set of 
parameters describing the system Hamiltonian which are 
under the observer's control. This is motivated by practi- 
cal considerations, since it is as of yet easiest experimen- 
tally to externally control aspects of the Hamiltonian. In 
this case the feedback step involves no measurement, and 
the observation and feedback processes may be regarded 
as qualitatively different, as in the classical theory. In 



view of this, the polar decomposition theorem motivates 
some definitions. 

By Kraus's representation theorem p3], every valid 
quantum evolution (a quantum operation) may be writ- 
ten as a POVM given by a set of operators fi n , where the 
probability of each outcome is P(n) = Tr[f2TO n p] an d the 
state resulting from each outcome is p n — Q, n pQ,' n /P(n). 
The only constraint on the f2 n 's is that ^] fit fl„ = 1. 
However, from the polar decomposition theorem, each 
of the operators Q may be written as the product of a 
unitary operator and a positive operator, so that 



(13) 



This provides a natural decomposition of a general quan- 
tum operation in terms of measurement and feedback. 
Consider first the action of the unitary operators. By 
themselves they do not describe the acquisition of infor- 
mation, and in that sense they do not describe a measure- 
ment. This can be seen from the fact that a unitary op- 
erator does not change the von Neumann entropy of any 
state it acts upon, and consequently extracts no informa- 
tion. However, unitary operations are precisely the kind 
that can be applied by Hamiltonian feedback. Hence, the 
unitary operators appearing in the polar decomposition 
may be thought of as characterizing purely the feedback 
part of the quantum operation. Note that we have writ- 
ten the polar decomposition so that the action of the 
unitary operator follows after the action of the positive 
operator, being a necessary condition for feedback. 

Conversely, the positive operators characterize the ac- 
quisition of information. They may always be written 
as a weighted sum of projectors, and therefore thought 
of as providing partial information about the states in 
the basis in which they are diagonal. When they corre- 
spond to rank 1 projectors, they provide complete infor- 
mation, in that the final state is pure. Since the unitary 
part has been factored out to obtain the positive opera- 
tors, we may regard these operators as representing pure 
measurement; the change induced in the quantum state 
is only that which is strictly necessary in order provide 
the information obtained during the measurement. We 
note that this decomposition of measurements into uni- 
tary and positive operators has been considered before 
in the context of measurements of the first and second 
kind ||. 

From this it is clear that every quantum evolution can 
be realized by a measurement in which the measurement 
operators are positive, followed by a feedback step in 
which the Hamiltonian is chosen to depend upon the mea- 
surement result. We see that the observation of a single 
observable, considered in section [n[ corresponds to the 
special case in which all the positive operators forming 
the POVM are mutually commuting. 

Under the above definitions damping processes, such 
as cavity decay and Brownian motion, are not considered 
pure measurements; they are viewed as equivalent to a 
fixed combination of measurement and feedback. Since 
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the object of feedback control is to limit the deviations of 
a system from a desired state (or more, generally, from a 
particular evolution - which means merely that the target 
state changes with time), feedback control is essentially 
a damping process (toward the target state). 

The polar decomposition theorem therefore fits snugly 
with the structure of Hamiltonian feedback, but it is nev- 
ertheless important to realize that this is not the only 
feedback process that may be considered in quantum me- 
chanics. First note that the product of two positive oper- 
ators need not be positive. Hence the evolution resulting 
from a sequence of pure measurements as defined above 
will in general be equivalent to a single pure measure- 
ment followed by some Hamiltonian evolution (ie. both 
measurement and Hamiltonian feedback). This is an il- 
lustration of the fact that quantum measurements involve 
'active' transformations of the states, as opposed to the 
'passive' measurements of classical physics plj . 

Consider now the full evolution of a system under 
Hamiltonian feedback control in a single infinitesimal 
time step dt, with initial state p. Since all dynamical 
processes commute to first order, one can treat even con- 
tinuous feedback control as alternating steps consisting 
of measurement and feedback. This is consistent with 
the general approach of this paper which is to consider 
the two steps separately. The system evolves under its 
own 'free' Hamiltonian, Hq, (which in many cases will 
be the desired evolution), and is affected by a source 
of environmental noise, which can be described by the 
non-selective evolution generated by a POVM. The mea- 
surement is also performed, and the feedback evolution 
applied. For a given measurement result n, we may write 
the full evolution as 

p n = e -i(^+W Pn (J2 W 3 pW})P n e^ +H ^ dt (14) 
j 

where the tilde indicates that we have not bothered 
to normalize the final state, and Wj are the operators 
describing the (undesirable) action of the environment. 
Since all the operators always commute to first order in 
dt, we have combined the free Hamiltonian with the feed- 
back Hamiltonian in the exponential. The task of feed- 
back control is to choose operators P n and H n such that 
the evolution is closest to the desired evolution. Before 
we consider this for Hamiltonian feedback, let us exam- 
ine what can be done in the absence of the conditional 
unitaries, using measurement alone, and the difference 
between the two kinds of feedback. 

By the definition above, using measurement alone one 
is restricted to POVM's in which all the measurement 
operators are positive, along with some overall unitary 
evolution independent of the measurement results. Now, 
to evaluate the efficacy of the control procedure, we must 
have a 'cost function' which measures how well we have 
achieved the control objective, as discussed above in the 
introduction. Since we have a desired 'target' state, 
a = |V , t)(V , t|j in mind at some final time (to be achieved 



following a single measurement, or a series of measure- 
ments), sensible cost functions will provide a measure of 
how close the final state, pt, is to the target state. A 
number of measures are possible, such as the inner prod- 
uct (Tr[p f cr]), the fidelity (Tr[ A /cr 1 / 2 p f0 -V2]) ) r the des- 
tinguishability ((l/2)Tr[/>f — <j]). Since we are interested 
only in target states which are pure, the Fidelity is simply 
the square root of the inner product, so that they pro- 
vide equivalent optimization problems. Throughout this 
paper we will use these as the quantities to be optimized. 

Now, the final state resulting from averaging the re- 
sults of a single pure measurement is given by 

Pi = J2 P nPPn- (15) 

n 

Since P n = P£, Ando's result (2fJ states that pi is al- 
ways majorized by p, which means that the eigenvalues 
of pi are at least as evenly distributed as the eigenvalues 
of p. This means that the von Neumann entropy of pf 
is always at least as large as the entropy of p. Another 
way of putting this is that each eigenvalue of pf is some 
weighted average of one or more of the eigenvalues of p. 

It follows almost immediately from the above results, 
that the fidelity of the final state cannot be any larger 
than the maximum eigenvalue, A max ( p ) , of the initial state 
p. To see this we first note that since all the eigenvalues of 
the final state, Xj are a weighted average of the eigenval- 
ues of p, none can be larger than the largest eigenvalue of 
p. Now, writing the fidelity in terms of the eigenvectors 
of P f, \4>j) we have 

(^tWW^^I^IWI 2 . (ie) 

3 

Since J^. |(<^j|^>t)| 2 = 1, the fidelity is merely a weighted 
average of the eigenvalues of pf, which proves the result. 
In fact, choosing any basis \ipi), we obtain the probability 
distribution over these states as 

Mi = {^i\pi\i>i) =^2xij^ji (I 7 ) 

3 

where Xij = K^ilV'i)! 2 - Since the matrix Xij satisfies 
YliXij = 1 an d 12 j Xij = 1, it is a doubly stochastic 
map, with the result that the vector {pti} is majorized 
by the vector {Aj}, and hence the von Neumann entropy 
of the distribution over any set of basis states is always 
at least as large as the distribution over the eigenvectors. 
Another way of saying this is that diagonal elements of 
a matrix resulting from a unitary transformation per- 
formed on a diagonal matrix are always at least as uni- 
formly distributed as the original elements (and almost 
always more so). 

Clearly this result for the upper bound on the final 
fidelity also holds for repeated measurements, in which 
subsequent measurements are not conditioned on the re- 
sults of previous measurements (i.e. for pure measure- 
ments with no feedback). However, it does not hold for 
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sequences of conditional measurements. In this case the 
initial state seen by subsequent measurements cannot be 
written as the state given by averaging over the results of 
previous measurements, since each final state may have 
a different measurement performed on it . 

It turns out that if we allow ourselves an infinite mea- 
surement strength, then the upper bound on the final en- 
tropy derived above can always be achieved in the limit 
of an infinite number of measurements. To see this one 
simply follows the procedure of Aharonov and Vardi, re- 
ferred to as the 'inverse quantum Zcno effect', developed 
in reference [p7f . Consider first an initial pure state, \ip), 
and the projector, P e = \e)(e\ onto the state 



(cos(e)\4>) + sin(e)|V*r)) 



(18) 



For e — this is the initial state, and for e — ir/2 this 
is the final state. For any value in between this state 
represents a rotation through an angle e from the initial 
state to the final state. If P e = \e)(e\ makes up one of the 
operators describing the measurement, the probability of 
failing to obtain the state \e) is 



P(e) 



He) 



(19) 



If we succeed in obtaining the state |e), then in that mea- 
surement step we have succeeded in rotating the state 
through e toward the desired state. We can attempt to 
rotate the state through the full n/2 radians by choos- 
ing e = 7r/(2M) and using M measurements. Since the 
probability of failure at each step is then second order in 
(1/M), while the number of steps scales only as M, as 
the number of steps tends to infinity, the total probabil- 
ity of failure tends to zero, and in this limit one achieves 
the desired rotation. To see that this achieves the upper 
bound when the initial state is mixed, we choose the pro- 
jector so as to rotate the eigenvector of p corresponding 
to the maximum eigenvalue to the desired state. For each 
measurement the corresponding POVM is then given, in 
general, by P ei & — 1 ; where P £i is the projector 
for the i l h measurement, and the f2; are arbitrary. The 
final state may then be written 



M 



M 



i i 

C e . [A] = P et AP £z + ttiAfl] 



N 



.n=2 



(20) 
(21) 



where p\n) = A n |n), with Ai the maximum eigenvalue, 
and A is an arbitrary operator. As the number of mea- 
surements tends to infinity, the first term on the RHS of 
Eq.(p0[) becomes 



M— >oo 



\\\ifa){ifa 



(22) 



Since this term contributes A max (p) to the fidelity, and 
since the other pure states making up the final density 



matrix cannot contribute negatively, the upper bound is 
achieved. 

What happens when we allow ourselves infinite mea- 
surement strength, and sequences of measurements in 
which subsequent measurements are conditioned on pre- 
vious results (i.e. measurement-only feedback)? In that 
case it is clear that the desired state can always be ob- 
tained with certainty; one begins by making a projection 
measurement in an arbitrary basis, which results in a 
set of pure states. Then the above procedure is used to 
rotate the resulting state to the target. 

When infinite strength is available, measurement-only 
feedback is equivalent to Hamiltonian feedback, since 
both allow any state to be created. However, in many 
continuous feedback control applications, the strength of 
the measurement is unlikely to be so much stronger than 
either the environmental noise or the free system dynam- 
ics that it can be used in this fashion in place of Hamil- 
tonian feedback @,||,|9). With measurements of finite 
strength the outcomes are necessarily random, so that 
Hamiltonian feedback cannot be simulated reliably. One 
can expect therefore that real applications will find the 
use of Hamiltonian feedback invaluable. 



V. MEASUREMENT: MAXIMAL INFORMATION 
AND MINIMAL DISTURBANCE 



In Section 111 we introduced the concept of the infor- 



mation provided about the system by a measurement, 
and this involved specifically the information regarding 
the state resulting from the measurement, a definition 
motivated by feedback control. This in turn motivated 
the definition of the strength of a measurement (e.g. s v or 
s p ), important because it constitutes a natural constraint 
when considering the optimization of control strategics. 
However, the actual information provided by a given 
measurement is not only a function of the measurement 
strength, but also the state of the system immediately 
prior to the measurement. As a result, once the available 
measurement strength is known, one can ask how to opti- 
mize the information provided by the measurement given 
the current state of the system. This defines the concept 
of a measurement returning maximal information (for a 
fixed measurement strength). 

In addition to providing information, quantum mea- 
surements can also introduce noise, a statement which 
we will now make precise. Consider first a classical sys- 
tem driven by noise. One can characterize the extent of 
the noise in some time interval by the increase in the en- 
tropy of the phase space probability distribution for the 
system state which is given by averaging over the noise 
realizations. This tells us how much we expect the noise 
to spread out the system in phase space in that time 
interval, and characterizes our uncertainty about the fu- 
ture state of the system resulting from the noise. Now 
consider a classical measurement. Since the initial state 



G 



is uncertain (or else we would not need to make the mea- 
surement), the result of the measurement is random, and 
as a result one's state of knowledge changes in a random 
fashion. However, this random change should not be con- 
sidered noise, since if one averages over all the possible 
measurement results (all the possible random changes) 
the probability distribution for the state of system re- 
mains unchanged. This is the sense in which classical 
measurement introduces no noise into the system. 

Now consider a quantum system driven by noise. The 
equivalent of the phase space distribution is the density 
matrix. In the same manner one can characterize noise 
by the resulting increase in the von Neumann entropy 
of the density matrix resulting from averaging over the 
possible noise realizations. One can therefore character- 
ize the noise introduced by a quantum measurement by 
calculating the increase in the von Neumann entropy (or 
alternatively the decrease in the purity) of the density 
matrix which results from averaging over the possible 
measurement results. While we saw above that in the 
classical case the measurement introduces no noise, this 
is not, in general true for quantum measurement. In 
terms of the von Neumann entropy, the excess noise in- 
troduced by a measurement is 

iV e v = S y (pt) - S y (p). (23) 

Defining it in terms of the purity we have 

7V C P-Tr[p 2 ]-Tr[p f 2 ]. (24) 

This makes precise the intended meaning of our initial 
statement that quantum measurements can introduce 
noise. Note that this has nothing to do with the Heisen- 
berg uncertainty principle — what concerns us here is 
the uncertainty of the future quantum state, and not the 
uncertainty of some set of observables for a given state. 
Recall that this is because the object of the control is the 
state of the system, and it is up to the observer to decide 
what the desired state is. Whether it be a minimal un- 
certainty state in the sense of the Heisenberg uncertainty 
principle is immaterial. 

Let us consider first the question of minimizing the dis- 
turbance due to the measurement. Recall that for a pure 
measurement, the evolution given by averaging over the 
measurement results is given by Eq. (|l5|), where all the 
P n are positive operators. Once again invoking Ando's 
result, we have that the von Neumann entropy of the final 
state is never decreased by the measurement. Measure- 
ments which minimize noise are therefore the measure- 
ments which leave the von Neumann entropy unchanged. 
These measurements are in this sense most like classical 
measurements. A set of measurements satisfying this cri- 
terion are those in which all the P n commute with the 
initial density matrix. In this case we have 

Pi = PnpPn = J2 P nP = P- ( 25 ) 



In the language of continuous measurements, since the 
operators P n are diagonal in the eigenbasis of the ob- 
servable, this means choosing to measure an observable 
which shares an eigenbasis with the density matrix. 

On a practical note, for continuous observation, mea- 
suring in the eigenbasis of the density matrix involves 
continuously changing the measured observable (note 
that such a process has been considered previously in 
the context of adaptive measurements ]30| ) . In many sit- 
uations this flexibility may be only partially available, or 
not at all. However, the above analysis indicates that for 
the purposes of noise minimization, one should choose 
the measured observable to be that in which the system 
is diagonal, or nearly diagonal, for the longest time dur- 
ing the period of control. In fact, this introduces the 
possibility that in certain cases it may be desirable to 
turn off measurement for periods in which the system 
occupies states which have large off-diagonal elements in 
the eigenbasis of the observable. Of course, the resulting 
noise reduction would have to be balanced against the 
accompanying loss of information. 

Maximizing the information for a fixed measurement 
strength is a much more difficult problem. Here we will 
examine a specific example for the continuous measure- 
ment of a two-state system. In the formulation of con- 
tinuous measurements that was discussed in Section [n|, 
we used measurement operators where each was a sum 
over an infinite number of projectors. For a two state 
system it is possible to obtain the same result (i.e. the 
same continuous measurement driven by Gaussian noise) 
by using a formulation with only two measurement oper- 
ators where each is the sum over only two projectors. To 
obtain a continuous measurement of a given observable 
the POVM is given by fi 2 , + Sl\ = 1 where the measure- 
ment operators are 

fi = ^IOXOI + a/T - k|1)(1|, (26) 
Ox = v^llXll + vT^ «|0)(0| (27) 

in which k = 1/2 + y/kdt, and |0) and |1) are the eigen- 
states of the observable. In each time step dt this pro- 
duces one of two results. The sum of these, in a time 
interval At = Ndt in which N results are obtained, is nat- 
urally governed by the binomial distribution. In the limit 
of large N (and infinitesimal At) this tends to a Gaus- 
sian, and one obtains the measurement record (Eq.(^)) 
and SME (Eq.©) given in Section ||, where the mea- 
sured observable Q = |0) (0| — 1 1) (1 1 . We can alternatively 
think of this measurement as a single-shot measurement, 
and in that case k can take any value between zero and 
one. Note that when n = or n = 1 the measurement 
is one of infinite strength. As k becomes closer to 1/2, 
the strength reduces, and for k = 1/2 the measurement 
provides no information. 

We can obtain measurements of all possible observ- 
ables by applying to the measurement operators an arbi- 
trary rotation over the Bloch sphere, given by the unitary 
transformation 
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U(9, 0)|O) = cos(6>/2)|0) + sin(0/2)|l) (28) 
U(9, 0)| 1) = cos(0/2)|l) - e - ** sin(0/2)|O). (29) 

Recall that this unitary transformation of the measure- 
ment operators preserves the measurement strength as 
defined in section III. Without loss of generality, we can 



choose the initial density matrix to be diagonal, and write 
it as p = p|0)(0| + (1 — One can then obtain an 

analytic expression for the final average purity which is 
given by 



early on in the control process, then it may prove desir- 
able to obtain information more quickly, at the expense 
of introducing extra noise, since the large uncertainty 
will be the major factor in reducing the effectiveness of 
the feedback. However, once the observer's knowledge is 
sufficiently sharp, it may prove more effective to reduce 
the nois e at the expense of some added uncertainty. In 
Section VII we will present simulations to show how this 
information trade-off affects the performance of feedback 
control in a two-state system. 



^ Tr[untut P un n W] 



E 



(30) 



Information and Disturbance vs. Measured Observable 

0.87 1 1 1 1 1 , i - 1 1 1 1 1 0.0 



This expression is fairly complex, and we will not need it 
here. (For a detailed analysis of this expression, includ- 
ing analytic expressions for general two-outcome mea- 
surements on two-state systems, the reader is referred 
to |32[ ). It is explicitly independent of <j>, as one would 
expect, since it is 9 alone which gives the angle (on the 
Bloch sphere) between the basis in which the density ma- 
trix is diagonal, and the basis of the measured observ- 
able. The final average purity is then explicitly depen- 
dent on the three parameters, p, k and 9, and we are 
concerned with maximizing this with respect to 9. When 
the measurement is non-trivial (k ^ 0.5), the strength 
of the measurement is finite (n =/= 0,1) and the initial 
state is impure (p G (0, 1)), one finds that the location 
of the maximum is independent of p and k, and occurs 
for 9 = 7r/2. This means that on average the maximum 
information is obtained about the final state (for fixed 
measurement strength) when the basis of the measured 
observable is maximally different from the eigenbasis of 
the density matrix (ie. if the density matrix is a mixture 
of u z eigenstates, then one should measure a x or a y ). 

We see then, that at least for a two-state system, the 
minimal disturbance is obtained when the measured ob- 
servable has the same eigenbasis as the density matrix, 
and the maximal information is obtained when the two 
bases are maximally different. Thus we obtain the result 
that, at least for a two-state system, there is a trade-off 
between information and disturbance in quantum feed- 
back control (in contrast to classical feedback control). 
This trade-off for finite strength quantum measurements 
is also of interest from a purely fundamental point of 
view, and this is explored in detail in reference [3^]. We 
plot both the excess noise introduced by the measure- 
ment, and the average final purity resulting from the 
measurement as a function of 9 in Fig. [lj. For a fixed 
measurement strength one therefore has the choice be- 
tween choosing a measurement to minimize the noise, and 
consequently obtain better control of the system (in that 
the system will fluctuate less around the desired value), 
or obtain a more accurate knowledge of the system at the 
expense of increased noise. Which is most desirable may 
well depend upon the current state of knowledge. For ex- 
ample, if the state of the system is poorly known, perhaps 




FIG. 1. 



0.4 0.5 0.6 

cos 2 (8/2) 

Here we plot the information obtained about the 
final state (solid line), characterized by the final average pu- 
rity, if, and the excess noise introduced by the measure- 
ments (dashed line), N£ , against the measured observable, 
parametrized by the angle 9. The parameters are p — 0.1 and 
k = 0.75 



VI. OPTIMAL HAMILTONIAN FEEDBACK 

Optimizing the Hamiltonian feedback for each feed- 
back step consists in finding the feedback Hamiltonian 
which maximizes the fidelity with the target state, given 
the initial state of the system, which is fairly straight- 
forward. Before we consider optimizing an infinitesimal 
Hamiltonian, which is necessary for continuous feedback, 
let us find the optimal unitary transformation. This re- 
sult would be useful if a finite time were available for the 
feedback step, and feedback strength was not an impor- 
tant constraint, so that a 'rotation' through any angle 
could be performed by the feedback Hamiltonian. 

Denoting the state of the system at the beginning of 
the feedback step as p, the fidelity with respect to the 
target state at the end of the feedback step is given by 



F(p f ,<7) = Tr[vV/W/2], 



(31) 



where a = |V't)(V , t| is the target state, and U is the uni- 
tary transformation constituting the feedback. We wish 
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to find U to maximize F(pf, a) 
arbitrary A and unitary V, 



Observing first that, for 



max \Tr[AV] \ = max |Tr[\/ A^AV'V] \ 



max 
v 



= Tr[VIU] 



(32) 



where we have used the polar decomposition theorem for 
A (A — V ArAV), and the Cj(A) are the eigenvalues of 
VaTK. Setting A = p^a 1 / 2 , this gives 

F(jH,a) = max\Tr[(UpU^) 1/2 a 1/2 V}\ 
= max |Tr[[/p 1/2 [/V /: V]| 

<^M'^M' 2 03) 



where the final inequality uses the result by von Neu- 
mann pn}] . In the last line Xj(p) and Aj(cr) are the 
eigenvalues of p and ex respectively, ordered such that 
the largest eigenvalue of p multiplies the largest eigen- 
value of a, the second largest the second largest, and so 
on down to the smallest eigenvalue of both states. Now 
we need merely realize that we can achieve the upper 
bound by choosing U so as to diagonalize p in the ba- 
sis of <7, re-ordering the basis states such that the largest 
eigenvalue of p is attached to the eigenstate of a with the 
largest eigenvalue. Writing the eigenstates of a as \<jj) 
(with eigenvalues ordered by size), and those of p as \pj) 
(similarly ordered) , then the explicit construction for the 
optimal U is 



U 



3 



(34) 



For continuous observation, each feedback step acts only 
for an infinitesimal time, dt. In this case the final state, 
Pi, is given by 



Pi 



i[H,p}At-\[H, [H,p)]{Atf 



(35) 



where H is the feedback Hamiltonian. The fidelity of the 
final state with the target state is therefore 



(V*rNV*r) = (V"tHV"t) 

-i(i> T \[H,p]\i> T )At 



±(ih\[H, [H,p]]\^ T )(At) 2 + 



(36) 



The first term is fixed, so to optimize the fidelity in 
the infinitesimal time step we should maximize the coef- 
ficient of At, being the dominant term. If the target state 
commutes with p, then this term vanishes for all so 
that we cannot choose a Hamiltonian that will cause an 
increase in the fidelity that is first order in time. If this 
situation occurs only for vanishingly small times, then it 



will make effectively no difference to the feedback perfor- 
mance. However, in those special situations in which the 
intrinsic dynamics preserves the commutivity of pit) and 
the target state, it can be important to choose a Hamil- 
tonian which maximizes the term which is second order 
in time. One should note, however, that if one has free- 
dom to choose the measurement basis, one can always 
choose a basis which disturbs the state so as to break the 
commutivity of p(t) with the target state, eliminating 
the need to consider the second order term in the time 
evolution. 

The maximization must be performed under a reason- 
able constraint on the eigenvalues of H (i.e. a constraint 
that captures the concept of a limitation on the strength 
of feedback). A number of suitable constraints are pos- 
sible, such as a restriction on the maximum eigenvalue 
of H, the sum of the norms of the eigenvalues, the sum 
of the squares of the eigenvalues, etc. Here we choose to 
use the last of these constraints, namely 



< 



(37) 



To maximize the coefficient of At in Eq.(|36|), we first 
note that it may be written as the operator inner prod- 
uct Tr[fL4], where 



A = i\i/j T )(v\ - i\v){ipr\ 
\v) = p\ih) 



(38) 



The maximum of the inner product, under the condition 
that the norms of the operators are constrained, occurs 
when the operators are aligned: H = cA, where c is in 
general a complex number, but real in this case to pre- 
serve the Hermiticity of H . With this inner product the 
norm of H is 



Tr[H 2 ]=J2UH) 2 



(39) 



Naturally, we take the maximum value allowed under the 
constraint, setting Tr[i/ 2 ] = p. This fixes the magnitude 
of the proportionality constant c, and results in the fol- 
lowing explicit construction for the Hamiltonian 



H = ix[\">h){ih\>p]- 



where 



with 
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a = (iPt\p \4>t) 
b = (V>tM^t) 



(40) 



(41) 



(42) 
(43) 



It now remains to maximize the coefficient of (dt) 2 
in Eq. ([36]) . Recall this is only required under the con- 
dition that the first term is zero, which implies that 







pIV't) = AtIV't)- In this case the expression for the coef- 
ficient may be written as 

(^t|#(p-At/Wt>- (44) 

Now, denoting the eigenvalues of p by A„ (ordered in de- 
creasing order), the eigenstates by \n), and denoting the 
target state as the eigenvector with n — M, the above 
expression becomes 

N 

(ifr\H P H\fa) = ]T(A„ - A M )|(n|ff|M)| 2 . (45) 

n=l 

Now, the constraint on the Hamiltonian may be written 
Tt[H 2 } = Tr[^2 \n)(n\H^2\m)(m\H] 

n m 

= J2\(n\H\m)\ 2 = p, (46) 

nm 

being a constraint on the sum of the square magnitudes of 
the elements of H . Only the subset |(n|if|M)| 2 of these, 
with n ^ M, contribute to the expression to be maxi- 
mized. To obtain the maximum value for the expression, 
we must therefore set all the elements of H which do not 
contribute to it to zero, this allowing the contributing 
elements to be as large as possible. The constraint then 
becomes 

AT 

£ \(n\H\M)f = p/2, (47) 

n^M 

where the factor of one half is enforced by the Her- 
miticity of H. The expression can now be seen as an 
average of the eigenvalues of p over the 'distribution' 
P(n) = | (n\H \M) | 2 , which is normalized to p by the con- 
straint. The maximum value is therefore achieved when 
all the weight of the distribution is placed on the term 
with the largest eigenvalue. The solution is therefore 

\(l\H\M)\ 2 = \(M\H\l)f = p/2, (48) 

with all other elements zero. The explicit construction 
for the optimal H being 

F = v ^72(|1)(V>t| + |^t)(1|). (49) 

Note that this assumes that the target state is orthogonal 
to |1), being the eigenvector with the largest eigenvalue. 
If the |1) is the target state, then there exists no Hamil- 
tonian evolution which will increase the fidelity, since the 
fidelity is the maximum it can be given the current purity 
of p. In that case we are free to set H — for that time 
step. 

It is worth noting that since the magnitude of the 
feedback Hamiltonian and the strength of the continu- 
ous observation are uniformly bounded, the evolution of 
the system is continuous. Given this, since the Feed- 
back Hamiltonian is a continuous function of the system 



state, it is intuitively clear that the feedback algorithm 
is well-defined (and continuous) in the continuum limit 
for almost all sample paths. 

We now have a feedback algorithm which can be used 
in conjunction with a measurement strategy for feedback 
control. In the next section we will implement such a 
strategy for the control of a two-state system. 

VII. FEEDBACK CONTROL OF A TWO-STATE 
SYSTEM 

In the previous sections we have considered the mea- 
surement and Hamiltonian feedback parts of the control 
problem separately. This resulted in a straightforward 
choice for a Hamiltonian feedback algorithm, but did not 
result in a clear choice for the measurement strategy. 
This was because we were able to identify in the mea- 
surement process a trade-off between information and 
disturbance. Because of this, the optimal measurement 
strategy for a given application is likely to depend upon 
the relative strengths available for the measurement and 
feedback. For example, if the feedback Hamiltonian is 
relatively strong, then it is likely that it will be able to ef- 
fectively counter the disturbance introduced by the mea- 
surement, and therefore the measured observable should 
be chosen to provide maximal information and the ex- 
pense of maximal disturbance. When this is not the case, 
the most desirable measurement is likely to be that which 
introduces less disturbance at the expense of providing 
reduced information. 

To examine the performance of a feedback control al- 
gorithm we must run the algorithm many times in order 
to obtain the average behavior. This is computationally 
very expensive, and so we use massively parallel super- 
computers which are ideal for this task. The results we 
present here are obtained by averaging one thousand re- 
alizations of the control algorithm. 

To provide a simple example of feedback control we 
consider a spin-half system precessing in a magnetic field 
aligned along the z-axis. In the absence of any noise, a 
spin aligned originally along the x-axis would rotate at 
a constant angular velocity around the z-axis, and we 
take this to be the desired (target) behavior. To provide 
the control problem, we subject the spin to noise which 
dephases it around the z-axis (this could arise from fluc- 
tuations in the magnetic field). The master equation de- 
scribing the free (but noisy) evolution of the spin is thus 
given by 

p = -ihu[a z ,p\ - (3[a z , [a z ,p\], (50) 

where u> is the precession frequency in the magnetic field 
and (3 is the strength of the dephasing noise. To imple- 
ment feedback control we allow the observer (who is also 
naturally the controller) to measure the spin along an ar- 
bitrary spin direction v(t), with measurement constant 
k, and apply a feedback Hamiltonian, Hfb(t), obtained 
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using the algorithm presented in the previous section. 
The full evolution of the controllers' state of knowledge, 
including the measurement and feedback, is therefore 

dp = -ift[a z + H !h {t),p}dt - P[a z , [a z ,p]]dt 
-fc[er v(t) , [<j v ( t ),p]]dt 

+V2k(cr v{t) p + pa v(t) ~ 2Ti[cr v(t) p}p)dW (51) 

We now simulate the dynamics resulting from the feed- 
back control loop for different values of 9, being the an- 
gle between the eigenbasis of the instantaneous system 
density matrix and the instantaneous measured observ- 
able, as discussed in Section |y|. For these simulations the 
strength of the magnetic field is such that Tluj — tt, so that 
the spin rotates once in a time interval t = 1. The noise 
strength is (3 = 0.4, and the parameters for the control 
loop are k = 2 and feedback strength /i = 10. We start 
the system in a pure state with the spin pointing along 
the x-direction, and evolve the controlled dynamics for 
a duration of t = 2 (the purity and fidelity settle down 
to their steady state behavior by approximately t = 0.8). 
Averaging the fidelity and purity over the full length of 
the run, for different values of 9 we obtain figure |2|. Ex- 
amining the dependence of the purity on 9, we find what 
we expect from the discussion in Section |V[ That is, the 
average purity of the system increases with 9, achieving 
a maximum at 9 = tt/2. This reflects the fact that, on 
average, measurements with a larger value of 9 extract 
information from the system at a faster rate. 

Feedback-Control Fidelity and Purity vs. Measured Observable 



a. 

T3 

a 0.91 - 



"O 0.9 - 



9 (radians) 

FIG. 2. The average purity (diamonds) and fidelity (cir- 
cles) of a feedback control algorithm for different measure- 
ment angles 6. The parameters are precession frequency 
hui — n, measurement constant k — 2, noise strength P = 0.4 
and feedback strength p — 10. 

The behavior of the fidelity, in this case, is similar to 
that of the purity. As 9 increases, the feedback is suffi- 
cient to ensure that even though we can expect the noise 
to increase with 9, the increase in purity has more of an 



effect on the fidelity than the noise. The result is that, 
with these resources, it is best to choose 9 = tt/2 (so as to 
measure in a basis maximally different from that which 
diagonalises the density matrix). However, from our pre- 
vious analysis of the trade-off between information and 
disturbance, we cannot always expect this to be the case. 



VIII. CONCLUSION 

In this paper we have considered the problem of con- 
trolling a quantum system in real time using feedback 
conditioned on information obtained by continuous ob- 
servation. The question of how to effect the best con- 
trol given the system dynamics (including environmental 
noise), and constraints on available resources is highly 
non-trivial. Here we considered a simplified problem in 
which we examine the measurement process and the re- 
sulting Hamiltonian feedback separately. Our purpose 
was both to examine what concepts are motivated by 
the feedback control problem, and to explore the ques- 
tion of optimization in this simplified problem. A concept 
which arises immediately in considering feedback control 
is the strength of a measurement. This strength quan- 
tifies the amount of information which a measurement 
provides. Previous definitions of the information pro- 
vided by quantum measurements have focussed on in- 
formation regarding the state prior to the measurement. 
Here we have argued that it is the information regard- 
ing the state resulting from the measurement which is 
relevant to quantum feedback control, and introduced a 
concept of measurement strength accordingly. 

Since measurements disturb quantum systems, it is im- 
portant to understand how this relates to feedback con- 
trol. We showed how it is possible to quantify the con- 
cept of the noise introduced by measurements in a way 
that is relevant to feedback control. One finds that while 
classical measurements do not introduce noise, quantum 
measurements do, in general, although it is possible, at 
least in principle, to make continuous quantum measure- 
ments that are noise free. 

Having arrived at precise concepts of information and 
disturbance, we examined the special case of continu- 
ous measurements performed on a two-state system, and 
found that maximization of information, and minimiza- 
tion of noise were mutually exclusive goals, implying 
the existence of an information-disturbance trade-off in 
quantum feedback control. This highlights the complex- 
ity of the control problem. 

We also considered the Hamiltonian feedback part of 
the control process. Defining the cost function as the fi- 
delity with a target state, and the feedback strength as 
the norm of the Hamiltonian, we were able to obtain the 
Hamiltonian generating the optimal instantaneous feed- 
back. 

Here we explicitly consider control realized by choos- 
ing dynamics conditional upon a measurement process. 
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Because of this one can refer to this technique as us- 
ing a classical controller, since it works by taking a 
classical process (the measurement record) and altering 
the system Hamiltonian accordingly, all of which can be 
achieved by a classical system. It is therefore worth not- 
ing that, so long as we are considering the dynamics of 
the controlled system alone to be the important quantity, 
this is equivalent to control which is realized by connect- 
ing the system, via an interaction Hamiltonian, to an- 
other quantum system, where this second system is large 
enough to be treated as a bath [jl2j. In general, using a 
second quantum system in this fashion may be referred 
to as using a quantum controller. When the quantum 
controller is finite dimensional and restricted in its dy- 
namical response time, one can expect the performance 
of classical and quantum controllers to be somewhat dif- 
ferent, and this is an interesting area for future work. 
The question of how to best design feedback strategies 



to control noisy quantum systems is a complex one. How- 
ever, the study of this problem will help us to understand 
better how quantum measurement may be exploited in 
the manipulation of quantum systems, and as quantum 
technology advances we can expect that this question will 
become increasingly important in practical applications. 
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